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ABSTRACT 



We report on VERITAS very-high-energy (VHE; E>100 GeV) observations 
of six blazars selected from the Fermi Large Area Telescope First Source Catalog 
(1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the 
VHE band, taking gamma-ray absorption by the extragalactic background light 
into account. This allowed the selection of six bright, hard-spectrum blazars 
that were good candidate TeV emitters. Spectroscopic redshift measurements 
were attempted with the Keck Telescope for the targets without Sloan Digital 
Sky Survey (SDSS) spectroscopic data. No VHE emission is detected during 
the observations of the six sources described here. Corresponding TeV upper 
limits are presented, along with contemporaneous Fermi observations and non- 
concurrent Swift UVOT and XRT data. The blazar broadband spectral energy 
distributions (SEDs) are assembled and modeled with a single-zone synchrotron 
self-Compton model. The SED built for each of the six blazars show a synchrotron 
peak bordering between the intermediate- and high-spectrum-peak classifications, 
with four of the six resulting in particle-dominated emission regions. 

Subject headings: active galactic nuclei: general — gamma rays: individual(RGB 
J0136+361, RGB J0316+090, RGB J0909+231, RGB J1058+564, RGB J1243+364, 
RX J1436.9+5639) 
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1. Introduction 

Blazars are active galactic nuclei (AGN) with a relativistic jet pointed close to the Earth 
line of sight. AGN are thought to be powered by accretion onto a supermassive black hole 
(~ 1O 9 M ) at the center of the host galaxy and are characterized by a double-peaked spectral 
energy distribution (SED) in the uF u representation. 

The lower energy peak of the broadband SED is attributed to synchrotron emission of 
highly relativistic electrons and positrons in the presence of a tangled magnetic field. In lep- 
tonic models, the higher energy peak is produced via the inverse-Compton (IC) up-scattering 
by the relativistic leptons of the synchrotron photons (synchrotron self-Compton, SSC) or 
a photon field external to the jet (external Co mpton, EC). More d e tails regarding leptoni c 
non-thermal e missi on of blazars can be found in lDermer et al.l (119921 ); iMaraschi et al.l ( 119921 ); 
Sikora et al.l (119941 ) and the references therein. Alternative models associ ate the higher- 



energy peak to interaction of relativistic protons with an ambient photon field ( Aharonian et al. 
~T lBednarekll993l : bar et al.fl997l iMannheim boooi iMiicke fc Protheroejjoool- "" 



2002: 



Bednarek 


1993; 


Dar et al.lll997 


Mannheim 


2000; 



20001 ). or a hybrid population comprised of both leptons and hadrons fadttchej 120071 ) . 



Pohl fc Schlickeiser 



The blazar population is divided into two subclasses: flat spectrum radio quasars (FS- 
RQs) and BL Lacs. FSRQs are, to first order, more distant, more luminous, and have 
stronger emission lines than BL Lacs. Historically, BL Lacs have been sub-classified based 
on their radio and X-ray flux ratios as low-, intermediate- or high-frequency-peaked BL Lac 
objects (LBL, IBL and HBLs, respectively) as described in iPadovani fc Giommil ( 119951 ) and 
Bondi et al.l (120011 ) . More recently, a classification based o n the location of the synchrotron 



peak in frequency space {v synch ) has been proposed by lAbdo et al. I (j2010al ). with low- 
spectrum-peaked (LSP) BL Lac objects having v sync h below 10 14 Hz, intermediate-spectrum- 
peaked (ISP) BL Lac objects peaking between 10 14 and 10 15 Hz and high-spectrum-peaked 
(HSP) BL Lac objects showing a peak above 10 15 Hz. HSP BL Lacs are the most common 
extragalactic object to be detected at very-high-energies (VHE; E > 100 GeV), compris- 



i ng 33 of 41 VHE B 



(lAcciari et al. l2009al 



Hoffman 



2010; 



b 



Lac objects detect ed as of June 2012. There have a lso been 4 ISPs 



20101; 



Albert et al. 



Oriel l2010af ) and 4 LSPs (IMazin et al.l 120091 ; iMariotti 1 12011 
l20Q7h detected since the advent of VHE gamma-ray astronomy 
in the late 1960s. In addition, three FSRQs, three radio galaxies, and two starbu rst galaxies 
complete the catalog of associated extragalactic objects ( Horan &; Wakely 2008 1 \ 



The potential scientific impact of increasing the catalog of VHE emitting extragalactic 
objects is substantial. A significant fraction of the power released from these objects is within 



J http : //tevcat . uchicago . edu/ 
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the VHE band. This makes the measurement of VHE blazar spectra an important component 
of the overall understanding of these objects. With a better sample of well-measured VHE 
blazar spectra available for study, a population-based investigation of gamma-ray production 
in these objects through broadband SED modeling will be possible, providing means to 
answer the long-standing question of whether VHE gamma-ray emission results from leptonic 
or hadronic processes in AGN jets. We can also apply the model-inferred properties of these 
new discoveries to understanding how the gamma-ray production may differ among blazar 



subclasses, exploring the apparent blazar sequence an d evolution of the AGN ( jFossati et al. 
19981 : iGhisellini fc Tavecchiol [20081 iMever et al.lboilh . 



VHE blazars can also be used to constrain the optical to near-infrared extragalactic 
background light (EBL) density and evolution, as well as the nature of cosmic rays. The 
EBL encodes the integrated history of structure formation and the evolution of stars and 
galaxies in the Universe. Understanding these characteristics requires detailed theoretical 
modeli ng of all the processes that contribute, e.g . structure formation and stellar evo 



lution (jDominguez et al. 



Franceschini et al. 



2008 



2011 



Gilmore et al. 



Hauser fc Dwek 



20C 1 



2009 



Stecker et al. 



Primack et al l 120051 : iFinke et al.l 12010 
20061 According to some cos 



mic ray models, e.g. lEssey fc Warren I ( 120121 ). interactions of cosmic rays along the blazar 
line-of-sight can produce relatively hard VHE gamma-ray spectra as compared with the 
high-energy gamma-ray spectra, depending on the distance to the blazar. 

The current catalog of extragalactic VHE objects largely contains relatively nearby 
AGN; only three have a confirmed redshift above z = 0.3. The highest redshift blaz ar de- 



tected at VHE thus far is the FSRQ 3C279, at a redshift of 0.536 ( Albert et al.ll2008bh . The 
proximity of these VHE blazars is partly a result of gamma-ray absorption by the EBL. VHE 
gamma rays that propagate through the intergalact ic medium are absorbed by low energy 
EBL photons via p air production, 7 + 7 — > e + + 



f Nikishov!ll962l ; iGould &: Shreder Ill967 



Stecker et al.lll992l ). The absorption process deforms the VHE gamma-ray spectra emitted 



limits on the local density of the far-infrared EBL ( 


Stecker et al. 


1993 


2005; 


Aharonian et al. 


2006; 


Primack et al. 


2011; 


Orr et al. 


2011 


. 



The search for new VHE blazars is complicated by the fact that many of these objects 
do not yet have known redshifts. BL Lac objects, by definition, display very weak or no 
optical emission or absorption lines used for spectroscopic redshift measurements. Due to 
the interaction of VHE photons with the EBL, non-detection of blazars with no redshift 
information can be either attributed to the object being too distant or having an intrinsically 
low luminosity in the VHE band. This makes VHE discovery observations of blazars with 
no measured redshift a risky venture, although it has proven successful in the past, as 
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Acciari et al. 


2009a. 


2010; 


Aliu et al. 


2011) 



The small field of view of imaging atmospheric Cherenkov telescopes (IACTs, < 5.0°) 
makes new source discovery from a large-scale sky survey difficult and therefore the hunt for 
VHE emitting objects has historically involved targeted observations of source candidates 
selected fro m surveys at lower frequencies, such as the selection of hard X-ray candidates 
presented in ICostamante fc Ghisellinil ( 120021 ). In this way, VHE blazar candidate selection 
has relied on experimen ts such a s EGRET onboard t he Compton G amma Ray Observator y 
dThompson et al.l Il993h . Swift f lGehrels et aD bopj ) and ROSAT jTurriziani et all l2007h . 
The launch of Fermi in June 2008 has enhanced VHE blazar discovery programs, leading 
to new blazar disc overies including RBS 0413, RXJ0648 7+1516, 1ES 0033+595 a nd 1RXS 
J101015.9-311909 rtAliu et all [20121 boill ; TMariotti IfioilL lAbramowski et alJboii l. 

The Fermi Large Area Telescope (LAT) observes the entire sky in the energy range 
from 20 MeV to > 300 GeV every three hours and has better sensitivity than its predecessor 
EGRET. Within the first 11 months of operation, the 1FGL catalog r eported the detection 
of 1451 sources at a significance greater than 5a ( Abdo et al. 2010bh. A ma jority of these 
sources are, or are expected to be, associated with AGN (lAbdo et al.ll2010a ). It is known 
that a large majority of these blazars has not yet been detected by VHE instruments, as can 
be seen by the current TeV catalog, which only contains 51 associated extragalactic objects. 
However, the proximity of the Fermi-LAT energy band to that of the IACTs makes the 
1FGL catalog a good place to search for candidate VHE blazars. We report on VERITAS 
observations above 100 GeV of six candidate blazars. For the first time, a multiwavelength 
description of their SED including radio, UV, X-ray and gamma-ray frequencies is assembled 
and modeled using a SSC model. 



2. Target Selection 

The energy coverage of IACTs overlaps with that of Fermi-LAT above 100 GeV and 
extends to tens of TeV. Sources that are most likely to be detected by both Fermi-LAT 
and IACTs have high fluxes and hard photon indices in the Fermi-LAT energy band. A 
selection process was established to identify the best VHE candidates within the first Fermi- 
LAT catalog. Sources at low Galactic latitude (|b| < 10°) were excluded, with the intent of 
removing the majority of Galactic sources from the selection. Additionally, sources with low 
integrated flux above 100 MeV (FiooMeV < 2 x 10~ 9 ph cm _2 s _1 ), with soft photon indices 
(T > 2.0), or with a low number of associated photons {N pre d < 20), were excluded from the 
selection. Just over 200 sources passed this initial set of cuts. 
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The remaining sources were re-analyzed using the same data set as that used for the 
1FGL catalog to search for the presence of curvature in their spectra. A log-parabola pa- 
rameterization for the spectra was chosen as the alternative to the power-law model (null 
hypothesis). In addition, the data were analyzed only using data above 1 GeV to confirm 
that the results from the power-law fit found for the entire Fermi-LAT energy band agreed 
well with the power-law fit above 1 GeV. Sources showing significant spectral curvature or 
softening of the spectrum above 1 GeV with an improvement in likelihood value correspond- 
ing to 3<r were excluded from the final selection as these sources are not expected to exhibit 
bright TeV flux levels. For sources with spectra best fit by a power law that also matched 
the power-law fit above 1 GeV, the Fermi-LAT spectrum was extrapolated up to the VHE 
regime ( 150 GeV - 1 TeV). When a reliable redshift measurement was available, the EBL 



model of iFranceschini et al.l ( 120081 ) was used to estimate the extrapolated VHE gamma-ray 
flux. Otherwise an assumed redshift of z = 0.2 was used, a conservative value considering 
that most known TeV sources have redshifts less than z = 0.2. Finally, the sources were 
ranked based on their extrapolated flux in the VHE regime. Six VHE candidates, all BL Lac 
objects, were selected for observations with the VERITAS telescope between September 2009 
and June 2010 for 10 hours each, corresponding to the VERITAS 3% Crab flux sensitivity 
exposure timescale. 

Table 1 shows the VHE extrapolated integral flux from the 1FGL catalog power-law 
fits used for selection of the six candidate VHE emitting blazars. These relatively high 
integral flux values above 150 GeV, shown in percentages of the Crab Nebula AuXlI are 
shown in comparison to the updated power-law fits from 29 months of Fermi-LAT data for 
the steady sources, and data from the VERITAS coincidental window for variable sources. 
The analysis of the extended Fermi-LAT dataset (2008 August 4 to 2011 January 4) is 
detailed in section 4.2. These updated extrapolations show much lower expected integral 
flux values for each of the six candidates, reflecting the fact that as more LAT data were 
collected, better high-energy statistics showed the candidates to be softer and/or dimmer 
than found with the elevent month data set used in the 1FGL. Additionally, we compare 
these extrapolated values to the upper limits derived from the VERITAS observations, where 
the analysis leading to these upper limits is detailed in Section 4.1. 



2 Flux calculated according to the curved power law presented in 



Albert et al. 



(l2008ah 
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3. The Targets 



RGB J0136+391 (1FGL J0136.5+3905) was discovered in the third Bologna sky sur- 
vey of 408 MHz radio objects (IFicarra et al.lll985l ). It was later detected as a X-ray bright 



active galaxy in the Northern ROSAT all-sky survey terinkmann et aL 



1997). and identified 



Laurent- Muehleisen et al 



spectr oscop ically as a BL Lac object, w i th an IBL sub-classification from 
( 119981 ) and lLaurent-Muehleisen et all (119991 ). respectively Based on the optical and X-ray 
spectr al properties, this blazar was proposed as a VHE candidate blazar by lCostamante fc Ghisellini 
( 120021 ) . under the assumption that the redshift was less than 0.2. Bright gamma-ray emis- 
sion above 1 GeV was detected from this sourc e by Fermi after t hree months of operation. 



Only an upper limit below 1 GeV was reported ( Abdo et al. 2009j ) 



This blazar also showed 

constant emission in the first eleven months of LAT operation (lAbdo et al. Il2010bl ). This 



blazar remains without a known redshift, with no previous spectroscopic redshift measure- 
ments found in the available literature. We attempted a spectroscopic redshift measurement 
on 2009 September 17 (MJD 55091) using the Keck Low Resolution Imaging Spectrometer 
(LRIS) instrument, but measured only a featureless power-law spectrum that is characteris- 
tic of BL Lac objects (see Figure 1 upper-l eft panel). Detailed spectroscopic analysis of this 
data can be found in Kaplan et al. I (120121 ) . 



RGB J0316+090 fl FGL J03 16 .1+0904 ) was first detected by the NRAO Green Bank 
91 met er radio telescope (IBecker et al.lll99ll ). It was late r optically identified as a BL La c 
object (IFischer et al .1119981 ) and sub- classified as an IBL by lLaurent-Muehleisen et al.l ( 119991 ). 
No spectroscopic redshift measurements had been made prior to this work. We attempted a 
spectroscopic redshift measurement on 2011 March 5 using the Keck Echelette Spectrograph 
and Imager (ESI) instrument, resulting in a featureless power-law spectrum except for two 
unidentified absorption features (see Figure 1 u pper- right panel). Th e detailed spectroscopic 



analysis of the Keck ESI data can be found in iKaplan et al. I (120121 ) . 



RGB J090 9+231 (1FGL J090 9.2+2310) was first detected by the NRAO Green Bank 
radio telescope feecker et al. 1991 ) and was later classified as a radio- loud active galaxy 



later ( Mickaelian et al 



(iBrinkmann et al.lll997 ). The BL Lac optical counterpart was identified nearly a decade 



20061 ). The redshift reported by the NASA Extra galactic Database 



(NED E could not be found within the corresponding reference (z = 0.231. iBrinkmann et al.l 



( 120001 )) and was therefore taken as unknown. Inspection of publicly available Sloan Digital 
Sky Survey (SDSS) data revealed two Mg II absorption lines in the otherwise featureless 
optical spectrum. Assuming that these lines could be intrinsic, or due to the absorption by 
an intervening cloud, a lower limit on the redshift of z > 0.43 is derived (see Figure 1 middle- 



3 http : //nedwww. ipac . caltech . edu/ 
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left panel). This is the only source out of the six selected which has a neighboring Fermi- 
LAT detected blazar within the VERITAS 3° field of view. 2FGL J0910.9+2246, associated 
with TXS 0907+230, is located 0.61° aw ay from RGB J0909+ 231. The high redshift blazar 
TXS 0907+230 (z = 2.66 according to lliealev et al~l toO^ )) shows no signal in the 10- 



100 Ge V band, with an up per limit of 9.6x10 photons cm s reported in the 2FGL 



catalog ( iNolan et al. 1120121 ). An exclusion region of radius 0.3° centered on the blazar was 
nevertheless used in the VERITAS analysis to avoid possible contamination. 

RGB J1058 +564 (1FGL J1058. 6+5628) was first detected in the 6C radio survey 
( IHales et al.lll990f ). It was identified as a BL Lac object during the a ssociation of t he RO SAT 
all-sky survey with the Hamburg Quasar Survey (INass et al.lll996l ). iBondi et al.l (]200lf ) clas- 
sified the object as an IBL based on the optical and X-ray chara cteristics. This bla zar was 



detected within the first t hree months o f Ferm i-LAT operation (lAbdo et al. 1120091 ) with a 



broadband SED shown in lAbdo et al. I (12010a ) which does not include any VHE informa- 
tion. The blazar also shows a high level of flux v ariability (probability of variability: 79%) 
in the high-energy gamma-ray band, as shown in lAbdo et al. I (j2010bl ) and is the only one 
out of the six targets in this paper that has a redshift measurement. The absorption lines 
corresponding to the redshift 0.143 can be seen in the SDSS spectrum shown in Figure 1 
(middle-right panel). 

RGB J 1 243+364 (1FGL J1243. 1+3627) was first reported in the B2 catalog of ra- 
dio sources flColla et al.l Il973[). This target w as also determined to be a radio-loud ac- 



tive galaxy by ROSAT (jBrinkmann et al.lll997l ). and specifically classified as a BL Lac in 
Appenzeller et al.l ( 119981 ) . NED cited the SDSS data for a spectroscopic redshift of z — 1.065. 
Inspection of this publicly available SDSS data revealed no lines suggesting this redshift but 
instead revealed Mg II absorption lines that translate to a lower limit of z > 0.485 (see Fig- 
ure 1 lower-left panel). A recent photometric redshift of z= 0.5lo 12 from iMeisner fc Romani 
( 120101 ) is in agreement with the SDSS lower limit. 



RXJ1436. 9+5639 (1FGL J1437.0+5 640) was detecte d in the ROSAT all-sky sur- 
vey and identified as a BL Lac object by INass et al.l ( 119961 ). This HBL, as classified by 
Nieppola et al.l (120061 ) based on the frequency of the synchrotron peak, remains without a 
redshift, although the redshift of the ga laxy cluster within the same region of the sky is 
known to be z= 0.15 ( jBauer et al.l 120001 ). Inspection of the publicly available SDSS data 
shows a featureless spectrum, shown in Figure 1. 
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4. Multiwavelength Observations and Analysis 



4.1. VERITAS 



The VERITAS observatory is an array of four 12-meter diameter IACTs, located in 
southern Arizona. VERITAS is sensitive to photons between 100 GeV and several tens of 
TeV with an energy resolution of better than 20%. The instrument has a 5a point source 
sensitivity of 1% of the Crab Nebula flux in less than 30 hours with an angular resolu t ion o f 
less than 0.1° for a Cra b-like source with a spectral index of 2.5. See IWeekes et al.l (120021 ) 
and iHolder et al.l ( 120061 ) for a detailed overview of the instrument. 



The VERITAS observations of the six VHE candidate blazars were completed between 
September 2009 and June 2010 (MJD 55122-55383). These observations were taken in wobble 
mode, with an offset of 0.5° from the source position in each of four cardina l direc tions 
to allow simultane ous background measurement, as explained in iFomin et al.l ( 119941 ) and 
Berge et al.l ( 120071 ). The radio location of the counterparts as specified by NED were used 
for source position. 

Air shower even ts initiated by gamm a and cosmic-rays are reconstructed following the 
procedure outlined in lAcciari et al.l (120081 ) . The recorded shower images are parameterized by 
their principal moments, giving an efficient method for suppression of the far more abundant 
cosmic-ray background. Any events with a total charge less than 50 photoelectrons are 
removed from the analysis. A set of cuts is then ap plied to the parameters in order to 
reject background events (see details of this method in iKrawczynski et al. I (120061 ) ) . These 
parameters reject cosmic-ray like events having the mean scaled width and length of the event 
camera image smaller than 1.1 or 1.4, respectively. Additionally, the reconstructed altitude 
of the maximum Cherenkov emission is required to be higher than 8 km above ground level. 
Gamma-ray like events are extracted from a signal region with a radius of 0.14°, centered at 
the coordinates of the candidate source. 

The quality-selected livetime collected for each target ranges from 4.4 to 14.2 hours 
and results in no detections , with significances ranging from -1.1 to 0.9er, calculated with 
Equation (17) of iLi fc Ma I ( 119831 ). The VERITAS observations and analysis results are 
detailed in Table 2. Integral upper limits at 99% confidence level are calculated using the 
Rolke et al.l (120051 ) method, assuming a photon index of T = 3.0 for the differential power-law 
spectrum dN/dE = N(E/E )- r . This index was assumed as a moderately softened index as 
compared to the Fermi-LAT index range (1.74—2.9), expected due to absorption by the EBL 
and possible intrinsic turnover. Additionally, this index value is representative of a typical 
TeV blazar. Without detection, the real index in the VHE band remains unknown and 
any integral upper limit derived is dependent on the index assumed. The differential upper 
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limit is quoted at the decorrelation energy, the energy where the calculated flux has minimal 
dependence on the index. Changing the spectral index by ±0.5 changes the differential upper 
limits by less than 10%. 



Cogan 


( 


2008 


) and 


Daniel 


( 


2008) 



?he upper limits range from 1 to 3% of the 
integral Crab Nebula flux above the threshold energy. The energy threshold for each obser- 
vation is defined as the energy at which the differential rate of reconstructed gamma rays 
from the postulated source reaches its maximum and is accurate to within the 20% energy 
resolution of the instrument, a value that is dependent on the observation zenith angle and 
sky brightness. 



4.2. Fermi LAT 



The Fermi-LAT is a pair-conversion telescope sensitive from 20 MeV to > 300 GeV, 
which operates in survey mode. Further details about the characteristics and performance of 



the LA T Instrument on the Fermi Gamma-ray Space Telescope can be found in lAtwood et al. 



(120091 ). Presented here is the analysis of the Fermi LAT data for the six candidates described 
in Section 3. Although the targets were selected based on only 11 months of data, more data 
was available after the completion of VERITAS observations and this larger data set is used 
for the modeling. More specifically, the LAT data from the time period of 2008 August 4 
to 2011 January 4 (MJD 54682.7-55565.0) were used for the modeling analysis. Except for 
variable sources, the analysis procedure was identical for each of the sources and proceeded 
as follows. 

For each candidate, events were extracted from a region of interest (ROI) of 10° radius 
centered on the target coordinates. Events from the "diffuse class" with zenith angle < 100° 
and energy between 300 MeV and 300 GeV were selected. Only data taken during periods 
when the rocking angle of the satellite was < 52° were used to reduce contamination from 
the Earth limb gamma rays, which are produced by cosmic rays interacting with the upper 
atmosphere. The significance and spectral parameters were calculated using an unbinned 



maximum-likelihoo d method implemented in the LAT Science Tool gtlikd_| (ICash 1 11979 



Mattox et al.lll996l ). A background model was constructed including nearby gamma-ray 
sources and diffuse emission. All sources within 12° of the central source in the second 



Fermi-LAT catalog (2FGL, iNolan et al. I (120 121 )) were included in the model. The spectra 
of known pulsars were modeled by a power laws with exponential cutoffs. As in the 2FGL 



4 ScienceTools v9r20p0 with the post-launch instrument response function (IRF) P6_Vll_DIFFUSE. 
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catalog, a log-parabola function was used for sources with significant spectral curvature. 
Otherwise, spectra were described as a power law. The spectral parameters of the sources 
in the ROI were left free during the fitting procedure. Sources outside the ROI, but within 
the 12° range had their spectral parameters fixed to the 2FGL catalog values. The Galactic 
diffuse emission and an isotropic component, which is the sum of the extragalactic diffuse 
gamma-ray emission and the residual charged particle background, were modeled using the 
recommended files§. 

The LAT data were first analyzed to calculate the time-averaged gamma-ray flux and 
spectral parameters of each candidate. A second analysis was then performed to study the 
impact of the Sun, a bright gamma-ray source, on the flux of candidates located near the 
plane of the ecliptic (RGB J0316+090 and RGB J0909+231). Removing time intervals when 
the Sun was in the ROI of each candidate had a negligible effect on the analysis results. 

Spectral points and a light curve were calculated for each candidate, and a temporal 
analysis was performed to search for flux variability. The timescale of this analysis was 
adjusted based on the specific candidate flux levels. The flux in each energy or time bin was 
determined with the spectral indices of all sources fixed to the best-fit value s over the full 



energy and time interval. For an energy or time bin with a test statistic (TS; see lMattox et al. 



( 119961 )) less than 9 or fewer than 3 predicted photons (N prec i), a 95% confidence level upper 



limit was calculated. 

The light curves were analyzed to search for flux variability with a likelihood method 
assuming a constant flux for the null hypothesis, following the same procedure as used 
in the 2FGL catalog. Only two sources showed significant evidence of flux variability: 
RGBJ0316+090 and RGB J1058+564. For these two sources, a refined analysis was done, 
selecting time periods contemporaneous with the VERITAS observation windows. The dura- 
tion of this contemporaneous period was chosen such that a significant detection (> 5cr) could 
be attained, resulting in slightly extended windows of MJD 55055-55145 for RGB J0316+090 
and MJD 55160-55185 for RGB J1058+564 with respect to the VERITAS observation win- 
dow. For these candidates, a butterfly corresponding to the la confidence interval was used 
to represent the spectral information (Figure 3). 

The derived spectral indices for the differential power laws obtained for the candidates 
are relatively hard, ranging from 1.74 to 2.09, as compared to the 2FGL average spectral 
index of 2.21 ± 0.01. The integral fluxes above 300 MeV range from 1.55 to 12.4xl0 _9 ph 



5 The files used were gll_iem_v02_P6_Vll_DIFFUSE . f it for the Galactic diffuse and 
isotropic_iem_v02_P6_Vll_DIFFUSE.txt for the isotropic diffuse component available at 



http : / /f ermi . gsf c . nasa . gov/ ssc/ data/p6vll/access/lat/BackgroundModels . html 
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cm 2 s 1 , indicating that these six blazars are bright in the high-energy band. The detailed 
Fermi-LAT results for each of these hard-spectrum, bright BL Lacs are summarized in Table 
3. 



4.3. Swift XRT 



The X-ray Telescope (XRT) onboard the Swift satellite ( IGehrels et al.ll2004l ) is a focusing 
X-ray telescope sensitive to photons with energy between 0.2 a nd 10 keV. The data used 
for the broadband SED modeling were analyzed as described in iBurrows et al.l (120051 ) with 
the HEASof t package Version 6.9 and XSPEQj Version 12.6.0. All data were taken in photon 
counting mode and pile up effects are accounted for when count rates exceeded 0.5 counts 
per second through the use of an annular source region, with a 1-2 pixel inner radius and a 20 
pixel outer radius. Each observation is binned and fit with an absorbed power law between 
0.3 and 10 keV, with t he neutral hydrogen d ensity taken from the Leiden/ Argentine/Bonn 
survey of Galactic HI ( jKalberla et al. 1 120051 ) . 



X-ray variability is commonplace for both VHE detected and non-detected blazars. If 
more than one exposure exists for an object and no variability is detected, the de-absorbed 
power-law fit of the combined data set is used to constrain the SED modeling. If variability 
is observed between multiple exposures, results from these separate exposures are shown 
independently on the SED plot, as is the case for RGB J0136+391 and RGB J1058+564. 
These blazars show flux variability factors of ~2 and 3 between exposures, respectively. 
The multiple exposures taken on RGB J0909+231 do not provide sufficient statistics for 
application of an absorbed power-law model and are therefore summed before fitting. Using 
Cash analysis ( iCash 1119791 ) did not improve the fitting of the single low-statistics spectra. 
The summed exposure fit result is shown on the SED. 

Each of the absorbed power-law fits applied to the XRT data resulted in photon indices 
greater than 2, with 2-10 keV integral flux levels between 0.3 and 21xl0 _12 ergs cm _2 s _1 . 
The index values suggest that the synchrotron component peaks below keV energies, char- 
acteristic of ISP blazars. The analysis results for each observation are summarized in Table 
4. 



http : //heasarc . nasa . gov/ docs/ sof tware/lheasof t/ xanadu/xspec/XspecManual . pdf 
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4.4. Swift UVOT 



The Swift-XKT observations were supplemented with simultaneous UV OT exposures 
taken in the V, B, U, UVW1, UVM2, and UVW2 bands rtPoole et al.ll2008h . The UVOT 
photometry is performed using the HEASof t program uvotsource. The circular source region 
has a 5" radius and the background region consists of several 15" radii cir cles of nearby 



empt y sky. The results are reddening corrected using E(B-V) c oefficients (ISchleee 



19981 ). The Galactic extinction coefficients are applied according to iFitzpatrick I ( 119991 ). The 



et al. 



uncertainty in the reddening E(B-V) is the largest source of error, especially in the UV bands 
for blazars that have a large value of E(B-V). If more than one exposure exists in a specific 
band for an object, the data from the observation closest to the VERITAS exposure are 
used, although no significant variability is seen across any band for any blazar. A summary 
of the UVOT analysis results is presented in Table 5. 



5. Broadband SSC Modeling 



Leptonic models for blazar jet emission attribute the higher energy peak in the SED 
to the inverse- Compton scattering of lower energy photons off a population of non-thermal, 
relativistic electrons. These same electrons are responsible for the lower-energy synchrotron 
emission that makes up the first peak. The target photon field involved in the Compton 
upscattering can either be the synchrotron photons themselves, as is the case in SSC models, 
or a photon field external to the jet in the case for EC models. 

The previously described multiwavelength data are matched with archival radio data 
collected from NED. These data are used to test a steady-state leptonic jet model for the 
broad-band continuum emission from the blazars. Although it has been fo und that ISP 



BL Lac s are sometimes better represented by external- Compton models, e.g. lAcciari et al. 



( 12009bl ). taking into account the lack of redshift information and lack of constraints from the 
broadband data we prefer not to apply an EC model, which would introduce additional free 
parameters, compared to the SSC model applied here. The model-predicted flux reflects the 
absorption of VHE gamma rays by the EBL acc o rding to the redshift information summarized 
in Table 1 using the model from lGilmore et al.l (120091 ). The level of TeV absorption resu lting 
from this model is consis tent with the absorption derived from the iFinke et al.l ( 120101 ) and 
Franceschini et al.l (120081 ) models. 



The SSC mode l appl ied to the multiwavelen gth data is th e equili brium version of the 
Bottcher &: Chiangl (120021 ) model, as described in lAcciari et al.l ( 12009bl ). In this model, the 
emission originates from a spherical blob of relativistic electrons with radius R. This blob is 
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moving down the jet with a Lorentz factor T, which corresponds to a speed of fire. The jet 
axis is aligned toward the line of sight with an angle o b s , which results in Doppler boosting 
with a Doppler factor D = (T[l — (3-p cos ^obs]) -1 . In order to minimize the number of free 
parameters, we assume that 8 b s = often referred to as the critical angle, for which 

r = d. 

Within the model, nonthermal electrons are injected and accelerated into a power-law 
distribution Q{^) = Qoj~ g between the low- and high-energy cut-offs, 7 min and ^ max - The 
radiation mechanisms considered lead to an equilibrium between particle injection, radiative 
cooling and particle escape. This particle escape is characterized with an escape efficiency 
factor 77, such that the time scale of escape t esc = rj R/c. This results in a particle distribution 
which streams along the jet with a power L e . Synchrotron emission results from the presence 
of a tangled magnetic field B, with a Poynting flux luminosity of Lb- 

The two parameters L e and Lb allow the calculation of the equipartition parameter 
ese = Lb/L c . This equipartition parameter is used as an estimator of the feasibility of 
the model, where models which result in e^ e ~ 1 are preferred. If the particle energy 
density greatly dominates over the magnetic field energy density, namely a particle dominated 
scenario, then the magnetic field cannot serve to collimate the jet. Following this design, 
acceptable parameters should result in at least partition conditions with Lb > L e . 

The broadband SED for each blazar can be seen in Figure 3, with the SSC model 
parameters for each representation summarized in Table 6. For each blazar, the archival 
radio data are taken as upper limits as these measurements are believed to contain a large 
amount of radiation produced in the radio lobes in addition to the synchrotron emission from 
the jet. The modeling for each of the six blazars shows synchrotron peak locations ~ 10 15 
Hz, characteristic of borderline ISP/HSP blazars. 

RGB J0136+391: This blazar is modeled using the lower of the two variable X-ray 
states for three different assumed redshifts (z = 0.2, 0.3 and 0.4). The variability timescale of 
three months suggested by the factor-of-two variability observed between XRT exposures is 
not short enough to provide a constraint on the size of the model emission region. Under the 
assumption that the gamma-ray emission remained steady during the Fermi and VERITAS 
observations, the hard LAT spectrum and low VERITAS upper limit derived for this blazar 
suggest a steepening of the gamma-ray spectrum at E > 100 GeV, which could be caused 
by the EBL absorption if a redshift z > 0.4 is assumed. Alternatively, this apparent break 
could originate from uncorrelated variability in the high and very-high energy bands. More 
specifically, a low flux state in the VHE band during VERITAS observations could provide 
a redshift-independent explanation of the apparent spectral softening. Figure 3 shows the 
model predictions, corrected by EBL absorption, for each of the redshift values. Only the 
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model at z = 0.4 is compatible with the VERITAS upper limit, which also results in a 
framework with balanced radiation and particle energy. 

RGB J0316+090: Due to the variability detected in the high-energy band, this blazar 
is modeled with Fermi-LAT data which spans the complete time period sampled by shorter 
VERITAS observations. The UVOT errors resulting from the E(B-V) reddening correction 
for this blazar are so large that all exposures except the V band are unconstraining to the 
SED modeling and therefore not shown. Application of the SSC model for an assumed 
redshift of z = 0.2 results in a particle dominated framework consistent with the VERITAS 
upper limit. 

RGB J0909+231: The SSC model parameters used to describe the broadband emis- 
sion of this blazar are determined for a redshift of z = 0.5, based on the lower limit derived 
from Mg II absorption lines found in the SDSS data. The Swift X-ray data from three 
exposures have been summed in order to provide sufficient statistics for absorbed power-law 
fitting. The model results in a particle-dominated scenario with a slight discrepancy between 
the model and the Fermi upper limit in the 1 — 3 GeV energy bin. Moving the Compton 
peak accomodates this upper limit only results in an even less favored, particle-dominated 
emission state. 

RGB J1058+564: Due to the variability observed in the LAT band, the SED modeling 
for this blazar is done with LAT data from the time period coincident with the VERITAS 
observation window and for the lower of two X-ray states observed. The variability timescale 
suggested by the variability factor of three observed between XRT exposures is not short 
enough to provide any constraint on the size of the emitting region. Both X-ray states are 
shown on the SED for reference. The SSC model shows agreement with the broadband data, 
although results in a disfavored particle dominated scenario. 

RGB J1243+364: The modeling for this blazar is completed for a redshift of z = 0.5, 
based on the new lower limit found from Mg II absorption lines in the public SDSS data. 
The model agrees with the broadband data and permits parameters at equipartition. 

RX J1436. 9+5639: The modeling for the broadband data of this blazar allows marginal 
agreement with the Fermi data when completed for the redshift of the spatially coincident 
galactic supercluster, more specifically z = 0.15. The resulting model parameters are far 
below equipartition, suggesting a particle dominated jet. 
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6. Discussion 

Six promising TeV blazar candidates were selected from the 1FGL catalog for observa- 
tions with VERITAS. These observations resulted in no VHE detections. Initial selection 
of these candidates from 1FGL power-law fit extrapolation suggested integral flux levels be- 
tween 1 and 12% of the Crab Nebula flux above 150 GeV after accounting for absorption 
by the EBL. With additional Fermi-LAT data and more information about the blazar red- 
shifts, the expected fluxes were updated to levels between 0.3% and 3.5% of the Crab Nebula 
flux above 150 GeV. The VERITAS exposure times were allocated based on the initial ex- 
trapolated values, resulting in only one upper limit below the updated VHE extrapolation, 
namely the upper limit for RGB JO 136+391. The non-detection of this blazar suggests spec- 
tral steepening of the high-energy spectrum that can either be explained by intrinsic spectral 
curvature, redshift-dependent EBL absorption or uncorrelated variability. 

Multiband observations are presented and allow for the construction and SSC modeling 
of the radio through TeV broadband SEDs. The model applied to these blazars is suffi- 
cient to represent the broadband data for each of the six, with model parame ters roughly 



comparable to those found for other VHE-detected ISPs using the same model ( lAbdo et al. 



201 ll ; lAcciari et al.ll2009bl ). The magnetic fields obtained in the modeling are generally low, 
resulting in disfavored, particle-dominated jets. This condition could be relaxed by including 
an external photon fie ld for inverse - Comp ton scattering, allowing so lutions closer to equipar- 



tition, as was done in lAbdo et all ( 120111 ) and lAcciari et all (l2009bf ). However, this scenario 



has not been explicitly tested. The quality of the data sampling of the inverse- Compton 
peak and the fact that the redshifts are not well determined for the majority of the objects 
do not provide sufficient constraints for a model with the extra degrees of freedom associ- 
ated with adding an external source of seed photons. Similar parameters were also found 
for the borderline ISP/HSP TeV detected blazar PKS 1424+240, with the exception of the 
spec tral index for the in jected electron distribution, which was found to be very soft (q = 



5.1, lAcciari et all (120101 )). while the typical values for this model fall between q =2.3 and 
2.7. 

Each of these blazars shows a synchrotron peak frequency characteristic of the ISP/HSP 
divide, namely v S ynch ~ 10 15 . Additionally, each of these BL Lacs exhibits a similar flux 
level within the high-energy gamma-ray band, showing comparable ~L S ynch and Ljc- This 
commonality of sub-type and peak luminosities is likely a bias introduced to the selection 
process through the sensitivity of the Fermi LAT instrument being greater in the 300 MeV 
- 100 GeV band as compared to the >100 GeV band. Nearby blazars releasing a large 
fraction of power above 100 GeV are known to be good candidates for detection in the VHE 
regime. This type of power emission is characteristic of HSP BL Lacs, the most commonly 
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VHE-detected type of blazar, but the least frequently detected BL Lac in the 1FGL and 
2FGL catalogs. 

Although the selection of these six blazars from high-energy Fermi-LAT data did not 
lead to any new TeV blazar discoveries, the exercise has been very illuminating. We 
are reminded that when selecting candidates for observation with TeV instruments, ac- 
curate redshift information is crucial. Without this information it is difficult to deci- 
pher if the lack of TeV emission is due to the intrinsic emission mechanism or the ab- 
sorption of gamma rays by the EBL. The redshift lower limits that were found in the 
SDSS for RGBJ0909+231 and RGBJ1243+364 indicate distances where future TeV de- 
tection from a deeper exposure is unlikely. Continued TeV observations of RGB J0136+361, 
RGBJ0316+090 and RX J1436. 9+5639, the three sources that remain without solid red- 
shift information, could provide insight into the distance to these sources, while further 
observations of RBG J1058+564 can be directly applied to further investigate the emission 
mechanism at work within the blazar jet. 
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Fig. 1. — Optical spectra for the six BL Lacs selected from the 1FGL catalog and observed 
with VERITAS. The black shows the object spectrum, while the red shows the instrumental 
noise. Only one BL Lac had a confirmed redshift upon selection (RGB J1058+564, z=0.143), 
confirmed with the SDSS spectrum shown in the middle-right panel. Redshift lower limits for 
RGBJ0909+231 (z > 0.4305; middle left) and RGB J1243+364 (z > 0.485; lower left) are 
found from Mg II absorption lines in the SDSS spectra. A featureless SDSS spectrum is found 
for RXJ1436. 9+5639. A redshift measurement attempt for RGB J0136+391 using the Keck 
LRIS instrument shows a featureless power- law spectrum (upper left). The normalized ESI 
spectrum of RGB J0316+090 exhibits three unidentified absorption features (upper right). 



The spectral analysis of the Keck LRIS and ESI spectra are detailed in lKaplan et al. I ( 120121 ). 
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Fig. 2. — Fermi-LAT light curves, with units of 10 9 ph cm 2 s 1 , are shown for the six can- 
didate VHE-emitting BL Lacs for 29 months of LAT data (MJD 54682-55565; 2008 August 
4-2011 January 4). The beginning and end of the VERITAS observations are denoted by ver- 
tical grey lines in each panel. The short VERITAS observation periods for RGB JO 136+364 
and RGB J0316+090 can be seen, representing only 7 days each. Upper limits at 95% confi- 
dence level are shown for bins resulting in TS of less than 9, denoted by downward pointing 
black triangles. Only RGB J0316+090 and RGB J1058+564 show any significant variability. 
For these two sources, LAT data only for the time periods within the window of VERITAS 
observations are used to constrain the modeling. 
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Fig. 3. — Non-contemporaneous broadba nd SEP data for ea c h BL Lac shown with corre- 
sponding SSC modeling using the model of lBottcher &: Chiang! ( 120021 ) . The modeling results, 
corrected for the EBL absorption, are shown with grey lines. The archival radio data points 
are taken from NED and used only as upper limits, accounting for the fact that much of the 
radio emission may result from diffuse synchrotron emission in the radio lobes of the jet. See 
the text for a more detailed description of the model parameters, with values listed in Table 
5. 



Table 1. Summary of the high-energy Fermi-L AT power-law extrapolation beyond 150 GeV resulting from 
power-law fits from the 11 months of data used to select the candidates as well as the extrapolation based on the 
extended data set spanning 29 months. These extrapolated flux values factor in the gamma-ray absorption resulting 
from interaction with the EBL and are reported in % Crab Nebula flux units above the same energy threshold in 
order to allow direct comparison with the upper limits derived from the VERITAS observations. 



Counterpart 


Original z 


1FGL 


1FGL Integral Flux* 


1FGL Extrapolated Flux 


Updated z 


Updated Extrapolated Flux* 


VERITAS UL 


Name 


Used for Selection 


Indext 


>300 MeV 


> 150 GeV 


Lower Limit 


[%Crab] 


[%Crab] 








[xl0~ 9 ph cm^s" 1 ] 


[% Crab] 








RGB J0136+391 


0.2 


1.73 


9.5 


4.1 


0.2 


3.5 


1.7 


RGB J0316+090 


0.2 


1.72 


19.6 


8.1 


0.2 


0.3 


2.0 


RGB J0909+231 


0.2 


1.46 


6.6 


12.1 


> 0.4305 


0.5 


1.5 


RGB J1058+564 


0.143 


1.97 


15.7 


1.5 


0.143 


1.1 


1.1 


RGB J1243+364 


0.2 


1.74 


5.8 


2.2 


> 0.485 


0.8 


3.1 


RXJ1436.9+5639 


0.15 


1.45 


3.6 


11.7 


0.15 


0.7 


2.4 



to 

t Taken from lAbdo et al. ] boiobh . Index and flux extrapolated without error for target selection. 



*Computcd with power-law fit from 29 months of data for steady sources and VERITAS coincident time window for variable sources (RGB J0316+090 and RGB J1058+564). 



Table 2. Summary of VERITAS observations and analysis results. The significances are computed fro m counts 



extr acted from source and background regions (ON and OFF, respectively) according to Equation 17 of iLi &: Ma 



(119831 ) . The VHE 99% confidence level integral upper limits (ULs) are used for the modeling and are derived with an 
assumed photon index of 3.0. The percentage Crab Nebula flux values are given above the corresponding energy 

thresholds for each observation. 



Counterpart z VERITAS Observation ON OFF Q t Significance Thresholdt Integral UL Dccorrclation* Differential UL 

Name Livctimc [hr] Window [MJD] Events Events a Energy Above Threshold Energy [GcV] at Dccorrelation Energy 

[GeV] [Xl0~ 8 ph m^ 2 s~ 1 ] (% Crab) [x!0~ 7 m _2 s _1 TeV _1 ] 



RGB J0136+391 




!).!) 


55122-55129 


1422 


11224 


0.1277 


-0.2 


16S 


4.9 (1.7%) 


260 


la; 


RGB J0316+090 




4.4 


55126-55133 


698 


5560 


0.1282 


-0.5 


168 


5.7 (2.0%) 


275 


1.6 


RGB J0909+231 


>0.4305 


14.2 


55150-55204 


2141 


14310 


0.1499 


-0.1 


165 


4.3 (1.5%) 


260 


1.3 


RGB J105S+564 


0.143 


9.8 


55160-55185 


1415 


10405 


0.1417 


-1.1 


220 


2.1 (1.1%) 


350 


4.7 


RGB J1243+364 


>0.485 


11.5 


55189-55291 


1617 


11133 


0.1389 


0.9 


150 


9.9 (3.1%) 


250 


12.il 


RXJ1436. 9+5639 


0.15* 


13.0 


55275-55383 


1563 


12314 


0.125 


0.6 


240 


4.1 (2.4%) 


:',7r, 


8.9 



^Defined as the energy at which the differential rate of reconstructed gamma rays from the postulated source reaches its maximum, accurate to within the 20% energy resolution of VERITAS. 
■'"Normalization factor between size of the ON and OFF regions. i 

^The energy at which the upper limit is minimally dependent on the index assumed for the source. |\3 

I 1 I ^ 00 

R.cdshift of spatially coincident galactic cluster, taken from Bauer et al.l 120001). 



Table 3. Summary of Fermi-LAT observations and analysis results. A variability study was completed using 29 
months of data for all sources. For steady sources, the spectral analysis is completed for 29 months of data. The two 
variable sources RGB J0316+090 and RGB J1058+564 show spectral analysis results for LAT data coincident with a 
slightly expanded window as compared to the VERITAS observation window so as to allow a 5er detection. 



Counterpart 


1FGL 


Variability Detected? 


Prob* 


MJD Fit 


TS 


Index 


Integral Flux > 300 MeV 


Name 


Name 






Window 




r 


[xl0~ 9 ph cm- 2 s- 1 ] 


RGB J0136+391 


J0136. 5+3905 


no 


0.39 


54682.7-55565.0 


1430 


1.78+0.04 


12.0+0.7 


RGB J0316+090 


J0316. 1+0904 


yes 


0.00 


55055.0-55145.0 


45.7 


2.09+0.26 


9.5+3.3 


RGB J0909+231 


J0909.2+2310 


no 


0.06 


54682.7-55565.0 


123 


1.68+0.13 


1.8+0.5 


RGB J1058+564 


J1058.6+5628 


yes 


0.00 


55160.0-55185.0 


37.2 


1.98+0.24 


12.4+3.8 


RGB J1243+364 


J1243. 1+3627 


no 


0.60 


54682.7-55565.0 


627 


1.76+0.06 


5.6+0.5 


RX J1436.9+5639 


J1437.0+5640 


no 


0.24 


54682.7-55565.0 


132 


1.74+0.13 


1.6+0.4 



*Prob is the probability of steady emission as measured from ATS per degree of freedom. 



Table 4. Summary of Swift-XRT observations and analysis results . Phot on counting mode data were fit with 
absorbed power laws using HI column densities from iKalberla et al. I ( 120051 ). Unabsorbed data were used for SED 
modeling. For multiple observations showing no significant variability, exposures were combined to improve statistics. 
If variability was detected between observations, both results are shown on the SED, although the model is only 
shown for the XRT observation falling closest to the window of VERITAS observations. 



Counterpart 


Observation 


MJD 


Exposure 


HI 


Flux (2-10 keV) 


Photon 


x 2 (DOF) 


Used in 


Name 


ID 




Time [ks] 


[xl0 20 cm- 2 ] 


[xl0~ 12 ergs cm _2 s _1 ] 


Index 




SED 


RGB J0136+391 


00039107001 


55134 


1.6 


6.17 


12.1+1.0 


2.24+0.05 


38.0 (40) 


yes 




00039107002 


55223 


2.9 


6.17 


21.4+1.0 


2.16+0.03 


107.2 (84) 


yes 


RGB J0316+090 


00038370001 


54899 


2.6 


12.6 


7.0+0.6 


2.13+0.06 


34.7 (35) 


yes 


RGB J0909+231* 


00040540001 


55338 


1.3 










no 




00040540002 


55338 


1.9 










no 




00040540003 


55342 


1.6 










no 




combined 






4.63 


0.3+0.06 


2.7+0.2 


5.1 (5) 


yes 


RGB J1058+564 


00038215001 


54852 


3.8 


0.78 


2.2+0.2 


2.60+0.05 


42.4 (39) 


yes 




00038453001 


54888 


1.0 


0.78 


6.3+0.6 


2.48+0.07 


25.9 (27) 


yes 


RGB J 1243+364 


00038445001 


54875 


2.1 


1.4 


8.2+0.5 


2.33+0.05 


50.4 (46) 


yes 


RXJ1436. 9+5639 


00038435001 


54918 


5.1 


1.55 


2.4+0.2 


2.31+0.05 


40.6 (31) 


no 




00038289001 


54918 


1.9 


1.55 


2.4+0.2 


2.27+0.05 


43.3 (31) 


no 




combined 






1.55 


2.4+0.2 


2.29+0.04 


72.5 (78) 


yes 



*Statistics extracted from each single observation are too low to fit an absorbed power law. The combination of the three exposures is used. 
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Table 5. Summary of Swift-XJVOT observations and analysis results. 



Target 


Observation 


Date 


Band 


Frequency 


VF V 




Used in 




ID 


[MJD] 




[Hz] 


[Jy Hz 




SED? 


RPR TD1 3£iJ_'iQ1 


000301 n7nm 
uuuoy iu i uui 


^^1 34 
D0104 


TTVW1 
U V Wl 


1 1 4 v 1 0l5 
1 . 14 X IU 


(^ i^4-0 
\ i . id it u. ioj 


v 1 0l2 
X IU 


y 




00039107001 


55134 


UVM2 


1.34 x 10 15 


H 44 + sol 


x 10 12 


y 


rpr T03ifi-i_0Q0 


0003&370001 

UUUOOO I UUU1 


04oyy 


v 

V 


r. r.r Y i pil4 
O.OO X 1U 


(n OA 4- 1 04~\ 

1 Z.Z4 ZU. 1.U4 1 


x 10 12 


y 




0003&370001 
UUUOOO / UUU-L 


04oyy 


u 


Q K7 v i nl4 
O.O / X 1U 


( 9 Q9 4- 3 9.7\ 
^z. yz zc o.o / j 


x 10 12 






0003°.370001 
UUUOOO ( UUU1 


o4oyy 


TTVW1 
U V VV 1 


1 1 4 v 1 Ol5 
1.14 X _LU 


/o AA _|_ 7 
^Z.44 It ( .D4J 


v 1 0^2 
X IU 


* 

n 




0003P.370001 
UUUOOO I UUU1 


o'loyy 


T TVW9 
U V vvz 


1 48 x 10^ 


f9 °. 4- 1 8 ^ 
^Z.O HI lO. O ^ 


x 10^ 2 




RPR Tnono_i_OQi 

txxjxti juyuy-i-zoi 


UUU4U04UUU1 


OOOOO 


V 


o.oo X IU 


/o or _i_ n -\-\\ 
^z.oo it U. 11J 


\/ i nil 

X IU 


y 




00040^40001 


^338 

OOOOO 


g 


Q3 y 1 14 
u.yo x iu 


1 1 . t i ZC U.ZO j 


x 10 11 


y 




00040^40001 
UUUIUO'iUUUl 


OuOOO 




S ^7 v 1 0I 4 

O.O / X iu 


(n of\ 4- fl 91 ^ 

1 Z.ZU HI U.Z1 J 


x 10^ 


y 




00040^40001 
UUU'±Uu i ±UUUl 


OOOOO 


U V Wl 


1 14 X 10^ 


^1. i O X U. lO J 


x 10^^ 


y 




00040^40001 
UUU4U04UUU1 


r rooo 
OOoOo 


TTVM9 

U V 1V1Z 


1 34 v 1 0l5 
1.04 X IU 


/ 1 f;o _|_ n 99'\ 
l^i .DO it u.zz j 


v i oil 

X IU 


y 




00040^40001 

uuu4U04uuui 


OOOOO 


U V vvz 


1 48 v 1 0^5 
1.4o X _LU 


/i on _i_ n 1 
i .ou it u. ly ) 


v 1 OH 
X IU 


y 




00040^40009 

uuu4U04uuuz 


OOOOO 


TTVM9 
U V 1V1Z 


1 34 v 1 0l5 
1.04 X IU 


/l Qfi 4_ n 99~\ 

ii. yu it u.zz j 


v 1 OH 
X IU 


n 




00040^40003 
UUU4U04UUUO 


r rooo 
OOOOO 


U V vvz 


1 4°. v 1 0^5 
1.40 X IU 


(~\ 74 4- O 90^ 

(i . / 4 it u.zu j 


v 1 OH 
X IU 


n 


prip Tinro i 
nXjD J 1Uoo+od4 


UUUoozloUUl 


04oOZ 


V 


O.oo X IU 


( 1.4z ± U.U4J 


w i n"12 
X IU 


y 




UuuoozloUul 


KA Q c:9 

o4ooz 


1 T 

u 


O.O / X IU 


14.40 it U.U4J 


v 1012 
X IU 


y 




00038215001 


54852 


UVW2 


1.48 x 10 15 


(1.49 ±0.05) 


x 10 12 


y 




00038453001 


54888 


B 


6.93 x 10 14 


(1.64 ±0.07) 


x 10 12 


n 




00038453001 


54888 


u 


8.57 x 10 14 


(1.70 ±0.06) 


x 10 12 


n 




00038453001 


54888 


UVW1 


1.14 x 10 15 


(1.66 ±0.06) 


x 10 12 


n 


RGB J1243+364 


00038445001 


54875 


V 


5.55 x 10 14 


(1.28 ±0.05) 


x 10 12 


y 




00038445001 


54875 


13 


6.93 x 10 14 


(1.30 ±0.04) 


x 10 12 


y 




00038445001 


54875 


u 


8.57 x 10 14 


(1.42 ±0.05) 


x 10 12 


y 




00038445001 


54875 


UVWl 


1.14 x 10 15 


(1.34 ±0.05) 


x 10 12 


y 




00038445001 


54875 


UVM2 


1.34 x 10 15 


(1.49 ±0.05) 


x 10 12 


y 




00038445001 


54875 


UVW2 


1.48 x 10 15 


(1.47 ±0.05) 


x 10 12 


y 


RXJ1436.9+5639 


00038435001 


54918 


V 


5.55 x 10 14 


(9.03 ±0.09) 


x 10 11 


y 




00038435001 


54918 


B 


6.93 x 10 14 


(1.14 ±0.07) 


x 10 12 


y 




00038435001 


54918 


U 


8.57 x 10 14 


(9.71 ±0.05) 


x 10 11 


y 




00038435001 


54918 


UVWl 


1.14 x 10 15 


(1.06 ±0.06) 


x 10 12 


y 




00038435001 


54918 


UVM2 


1.34 x 10 15 


(1.21 ±0.07) 


x 10 12 


y 




00038435001 


54918 


UVW2 


1.48 x 10 15 


(1.12 ±0.06) 


x 10 12 


y 




00038289001 


54918 


V 


5.55 x 10 14 


(8.58 ±0.05) 


x 10 11 


n 




00038289001 


54918 


U 


8.57 x 10 14 


(1.06 ±0.04) 


x 10 12 


n 




00038289001 


54918 


UVW2 


1.48 x 10 15 


(1.26 ±0.06) 


x 10 12 


n 



*Galactic reddening dominates the uncertainty with E(B-V) = 0.356 for RGB J0316±090; U, UVWl, 
UVW2 band flux measurements do not constrain the SED. 



Table 6. Summary of the broadband SED SSC modeling parameters. See text for parameter descriptions. 



Counterpart Name 


z 


■v ■ a 

J in in 


imax 


</ a 


tesc [hr] a 


B [Gauss] a 


1 


Radius [cm] a 


L e [ergs/sec] a 


L B [ergs/sec] b 


Lb/L b 


b 




RGB J0136+391 


0.2* 


4.0 x 


10 4 


6.0 x 10 5 


2.6 


120 


0.15 


15 


3.4 X 10 16 


1.1 X 10 44 


2.3 x 10 43 


2.1 x 10 


-1 


25.9 


RGB J0136+391 


0.3* 


2.5 x 


10 4 


5.0 x 10 5 


2.6 


50 


0.4 


15 


2.4 X 10 16 


1.0 x 10 44 


7.8 x 10 43 


7.5 x 10 


-1 


19.2 


RGB J0136+391 


0.4* 


2.4 x 


10 4 


5.0 x 10 5 


2.6 


30 


0.45 


15 


3.0 x 10 16 


1.3 x 10 44 


1.5 x 10 44 


1.1 




25.9 


RGB J0316+090 


0.2* 


1.5 x 


10 4 


7.0 x 10 5 


2.5 


15 


0.11 


15 


8.0 x 10 16 


9.8 x 10 43 


6.5 x 10 43 


6.6 x 10 


-1 


59.2 


RGB J0909+231 


0.5 


2.3 x 


10 4 


5.0 x 10 5 


2.8 


10 


0.1 


20 


2.5 x 10 1S 


1.5 x 10 44 


9.4 x 10 42 


6.4 x 10 


-2 


17.4 


RGB J1058+564 


0.1433 


4.0 x 


10 4 


3.0 x 10 5 


2.9 


10 


0.25 


20 


9.0 x 10 15 


3.9 x 10 43 


7.6 x 10 42 


2.0 x 10 


-1 


4.8 


RGB J1243+364 


0.5 


3.0 x 


10 4 


1.0 x 10 6 


2.5 


30 


0.1 


20 


1.5 x 10 17 


3.2 x 10 44 


3.4 x 10 44 


1.1 




104.1 


RX J1436.9+5639 


0.15 


4.0 x 


10 4 


2.5 x 10 6 


2.6 


100 


0.01 


20 


2.7 x 10 17 


1.5 x 10 44 


1.1 x 10 43 


7.5 x 10 


-2 


143.7 



* Assumed rcdshift values due to lack of any redshift constraining measurement. 
a Parameters adjusted to match the model to the data. 
b Quantities derived from the modeling result. 



